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Preface 


The  numerical  water  quality  model  smdy  for  the  Los  Angeles  Harbor  Pier 
400  Project  was  conducted  at  the  U.S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  for  the  U.S.  Army  Engineer  District,  Los  Angeles. 

Mr.  Dennis  G.  Markle,  Chief,  Wave  Processes  Branch,  Wave  Dynamics  Divi¬ 
sion,  Coastal  Engineering  Research  Center  (CERC),  was  the  WES  study 
manager. 

The  study  was  conducted  by  Mr.  Ross  W.  Hall  of  the  Water  Quality  and 
Contaminant  Modeling  Branch  (WQCMB),  Environmental  Processes  and 
Effects  Division  (EPED),  Environmental  Laboratory  (EL),  WES.  Dr.  Harry 
Wang,  CERC,  coordinated  and  provided  the  hydrodynamic  model  output  data 
used  for  the  water  quality  simulations.  Dr.  Mark  S.  Dortch,  Chief,  WQCMB, 
provided  direct  supervision  and  report  review.  General  supervision  was  pro¬ 
vided  by  Mr.  Donald  L.  Robey,  Chief,  EPED,  and  Dr.  John  Keeley,  Director, 
EL. 

At  the  time  of  publication  of  this  report,  Dr.  Robert  W.  Whalin  was  Direc¬ 
tor  of  WES.  COL  Bruce  K.  Howard,  EN,  was  Commander. 

This  report  should  be  cited  as  follows: 

Hall,  Ross  W.  (1995).  “Numerical  Water  Quality  Model  Study 
for  the  Los  Angeles  Hartwr  Pier  400  Project,”  Miscellaneous 
Paper  EL-95-1,  U.S.  Army  Engineer  Waterways  Experiment 
Station,  Vicksburg,  MS. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


1  Introduction 


The  Port  of  Los  Angeles  plans  to  construct  an  additional  port  facility, 
referred  to  as  Pier  400.  The  Pier  400  harbor  facility  may  affect  water  quality 
by  changing  the  tidal  circulation  and  flushing  patterns.  The  objective  of  this 
study  was  to  compare  flushing  and  DO  resources  at  existing  conditions  and  at 
two  stages  of  plan  implementation.  The  comparisons  were  conducted  through 
numerical  model  simulations. 

The  numerical  water  quahty  model  (WQM)  used  for  the  comparisons  was 
developed  for  the  Los  Angeles  and  Long  Beach  Harbors  Model  Enhancement 
Program  (Hall  1990).  The  modeling  technology  consisted  of  three- 
dimensional,  time-varying  hydrodynamic  and  water  quality  models.  The 
hydrodynamic  model  (HM)  provides  the  circulation  required  by  the  transport 
terms  of  the  WQM.  The  WQM  is  indirectly  linked  to  the  HM.  Thus,  the  HM 
is  applied,  and  the  output  is  stored  and  subsequently  used  by  the  WQM. 

This  model  application  and  the  previous  model  enhancement  study  used 
water  quality  data,  meteorological  conditions,  and  tides  measured  during 
August  1987.  However,  the  model  enhancement  study  used  three  vertical 
layers,  whereas  this  model  application  has  five  layers.  Therefore,  to  demon¬ 
strate  that  the  model  could  be  used  for  the  Pier  4(X)  comparisons  using  five 
layers,  the  August  1987  conditions  were  included  in  the  tests.  The  differences 
between  the  1987  model  enhancement  study  condition,  the  existing,  and  the 
two  stages  of  plan  implementation  are  the  basin  morphologies.  Changes  in  the 
basin  morphology  result  in  changes  in  circulation  used  by  the  WQM.  In  this 
report  the  1987  condition  is  referenced  as  “Base,”  the  existing  condition  as 
“Existing,”  and  the  two  stages  of  Pier  400  implementation  as  “POLA  Stage  1” 
and  “POLA  Stage  2.” 

A  discussion  of  the  WQM  used  in  this  study  is  presented  in  Chapter  2. 
Chapter  3  describes  the  water  quality  and  flushing  comparisons,  and  Chapter  4 
contains  a  summary  of  results  and  conclusions. 
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2  Water  Quality  Model 


The  WQM  is  an  integrated  compartment  box  model  with  enhanced  horizon¬ 
tal  and  vertical  transport  schemes.  Interface  code  is  linked  to  the  HM  to  write 
time-invariant  data  and  temporally  averaged  flows.  The  water  quality  kinetic 
algorithms  were  adapted  from  the  HydroQual,  Inc.,  Potomac  Estuary  Eutrophi¬ 
cation  Model  (Thomann  and  Fitzpatrick  1982). 


Advective  and  Diffusive  Transport  Schemes 

A  modified  version  of  QUICKEST  (Quadratic  Upstream  Interpolation  for 
Convective  Kinematics  with  Estimated  Streaming  Terms)  scheme  (Leonard 
1979,  Hall  and  Chapman  1985)  is  implemented  in  the  WQM.  The 
QUICKEST  implementation  eliminated  excessive  numerical  diffusion  charac¬ 
teristic  of  upwind  schemes  or  numerical  instabilities  characteristic  of  centered 
difference  schemes. 

The  vertically  stretched  sigma  coordinates  in  the  hydrodynamic  model 
account  for  much  of  the  vertical  advection  through  layer  thickness  expansion 
and  contraction.  As  a  consequence,  the  vertical  velocities  are  small.  However, 
the  time  step  can  be  limited  due  to  vertical  diffusion.  Therefore,  an  implicit 
vertical  advective  and  diffusive  transport  scheme  is  implemented  in  order  to 
use  a  larger  time  step.  Central  differences  were  used  for  both  the  advective 
and  diffusive  terms. 


Model  Linkage 

The  HM  used  for  this  study  was  CH3D  (Curvilinear  Hydrodynamics  in 
Three  Dimensions).  The  hydrodynamic  and  water  quality  models  were  linked 
by  spatially  and  temporally  averaging  CH3D  output  to  drive  the  WQM.  The 
hydrodynamic  model  used  extensive  spatial  resolution  to  resolve  the  flow  field 
and  to  minimize  the  need  for  parameterization.  The  CH3P  spatial  resolution 
was  on  the  order  of  100  m  and  required  a  time  step  of  60  sec  for  stability.  In 
contrast,  the  WQM  has  characteristic  time  scales  on  the  order  of  hours,  which 
are  determined  by  the  kinetic  rate  coefficients.  The  desired  analyses  of  water 
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quality  allow  a  spatial  resolution  an  order  of  magnitude  less  than  that  used  by 
CH3D.  The  larger  horizontal  WQM  grid  size  allowed  the  use  of  a  time  step  of 
900  sec.  Both  CH3D  and  the  WQM  used  five  vertical  layers. 

Vertical  diffusivities  were  calculated  in  CH3D.  The  linked  interface  code 
temporally  averaged  the  calculated  values  and  included  the  averaged  vertical 
diffusivities  in  the  hydrodynamic  information  written  for  subsequent  use  by  the 
WQM.  However,  since  temperature  and  salinity  simulations  were  not  an 
objective  of  the  HM  task,  precise  tuning  of  CH3D  temperature  and  salinity 
calculations  to  provide  vertical  diffusivities  was  not  conducted.  As  a  result  the 
temporally  averaged  vertical  diffusivities  were  generally  0.001  m^  sec  *.  Verti¬ 
cal  diffusivities  of  0.001  m^  sec'*  in  the  WQM  resulted  in  uniform  vertical 
distributions  of  DO  in  the  WQM.  A  systematic  decrease  in  the  magnitude  of 
vertically  diffusivities  in  the  WQM  indicated  that  the  use  of  a  spatially  con¬ 
stant  value  of  0.0001  m^  sec'*  resulted  in  simulated  vertical  distributions 
matching  measured  values.  Field  dye  studies  conducted  by  Tekmarine  (1987) 
provided  an  order  of  magnitude  estimate  for  the  horizontal  diffusion  coefficient 
of  1  m^  sec  *  (Fisher  et  al.  1979). 

The  HM  grids,  WQM  grids,  and  WQM  grids  overlaying  the  HM  grids  are 
displayed  as  follows: 


Plate 

Plot 

1 

Base  HM  Grid 

2 

Base  WQM  Grid 

3 

Base  WQM  Grid  Overlaid 

4 

Existing  HM  Grid 

5 

Existing  WQM  Grid 

6 

Existing  WQM  Grid  Overlaid 

7 

POLA  Stage  1  HM  Grid 

8 

POLA  Stage  1  WQM  Grid 

9 

POLA  Stage  1  WQM  Grid  Overlaid 

10 

POLA  Stage  2  HM  Grid 

11 

POLA  Stage  2  WQM  Grid 

12 

POLA  Stage  2  WQM  Grid  Overlaid 

Kinetic  Routines 

The  WQM  kinetic  code  included  eight  variables  for  which  measured  data 
were  available  and  were  used  to  monitor  model  results.  These  variables  are 
referred  to  as  state  variables  and  include 
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a.  Temperature. 

b.  Tracer. 

c.  Phytoplankton. 

d.  PO4-P  (dissolved  phosphate). 

e.  NH4-N  (ammonium). 

/.  NO2+NO3-N  (nitrite  plus  nitrate). 

g.  CBOD5  (5-day  biochemical  oxygen  demand). 

h.  DO  (dissolved  oxygen). 

Temperature  is  a  primary  determinant  of  the  rate  of  biochemical  reactions. 
Tracer  was  a  conservative  state  variable  used  for  the  flushing  analyses.  The 
remaining  state  variables,  phytoplankton  through  DO,  were  transported,  inter¬ 
acted  with  each  other,  and  added  or  deleted  through  boundaries  or  sources  and 
sinks.  Phytoplankton  biomass  was  modeled  as  carbon.  However,  phytoplank¬ 
ton  biomass  was  reported  in  chlorophyll-a  (Tekmarine  1987).  A  carbon  to 
chlorophyll  (C/Chl)  ratio  of  80.0  was  used  to  convert  chlorophyll  to  carbon. 
Ph3^oplankton  are  referenced  as  algae  in  the  tables.  The  kinetic  routines  are 
described  in  HaU  (1990).  The  kinetic  constants  used  are  listed  in  Table  1. 

Salinity  influences  the  saturation  concentration  of  DO.  Salinity  was  set 
constant  at  32  parts  per  thousand  C/oJ-  A  constant  reaeration  coefficient  of 
1.0  m  day'*  was  used  for  the  water  quality  simulations.  The  application  of 
O’Connor’s  (1983)  relationship  between  wind  speed  and  reaeration  resulted  in 
a  daily  average  reaeration  coefficient  for  San  Pedro  Bay  that  varied  between 
0.96  and  1.02  m  day*  (Hall  1990). 


Initial  Conditions 


The  WQM  simulation  extended  from  August  1-31,  1987.  Initial  conditions 
were  specified  for  each  of  the  WQM  state  variables  and  for  each  cell  at  the 
initiation  of  the  simulation.  The  water  quality  data  collected  by  Tekmarine 
(1987)  on  August  4,  1987,  were  used  to  specify  initial  conditions. 

Examination  of  the  data  revealed  minor  horizontal  spatial  variability  and 
similar  vertical  distributions  among  stations.  Temperature  and  DO  were  mea¬ 
sured  at  2-in  intervals.  Chlorophyll-a,  CBOD5,  and  the  nutrients  PO4-P, 
NH4-N,  and  NO2+NO3-N  were  measured  1  m  below  the  surface,  middepth,  and 
1  m  above  the  bottom.  No  detectable  chlorophyll-a  or  CBOD5  was  measured 
August  4,  1987,  interior  to  the  breakwaters. 
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Table  1 

Kinetic  Constants 

Description 

Unit 

Solar  radiation  attenuation  constant 

0.50  m-’ 

Light  intensity  at  maximum  algal  photosynthesis 

70.0  W  m'^ 

Fraction  of  daylength 

0.5 

Half-saturation  constant  for  nitrogen 

0.025  g  N  m'® 

Half-saturation  constant  for  phosphorus 

0.001  g  P  m  ® 

Saturated  algal  growth  rate  at  20  “C 

2.0  day"' 

Temperature  coefficient  for  algal  growth 

1.068 

Algal  respiration  rate  at  20  °C 

0.1875  day^ 

Temperature  coefficient  for  algal  respiration 

1.045 

Algal  maximum  excretion  rate 

0.03  day-'' 

Algal  mortality  rate  at  20  "C 

0.02  day"' 

Algal  settling  velocity 

0.1  m  day"' 

Phosphorus  to  carbon  ratio 

0.025  g  P  (g  C)'' 

CBOD  deoxygenation  rate  at  20  "C 

0.2  day-'' 

Temperature  coefficient  for  CBOD  deoxygenation 

1.047 

Half-saturation  constant  for  CBOD 

0.5  g  Oj  m'® 

Ratio  of  ultimate  CBOD  to  CBOD5 

1.85 

Oxygen  to  carbon  ratio 

2.67  g  Og  (g  C)'' 

Nitrogen  to  carbon  ratio 

0.25  g  N  (g  C)’' 

Nitrification  rate  at  20  °C 

0.1  day'^ 

Temperature  coefficient  for  nitrification 

1.08 

Half-saturation  constant  for  nitrification 

2.0  g  Og  m"® 

Temperature  coefficient  for  SOD 

1.08 

Half-saturation  constant  for  SOD 

0.5  g  Og  m"® 

Because  of  the  observed  spatial  distribution,  sample  depths,  and  absence  of 
detectable  measurements,  the  following  procedures  were  used  to  assign  initial 
conditions: 

a.  Nondetectable  measurements  were  assigned  a  value  equal  to  one-half 
the  detection  limit  (Table  2). 

b.  Initial  conditions  were  assumed  horizontally  constant  but  varied  with 
depth.  Temperature  and  DO  measurements  were  assigned  to  a  layer 
based  on  the  depth  of  measurement.  For  the  constituents  measured  at 
three  levels,  the  surface  value  was  assigned  to  layer  1,  the  bottom  value 
assigned  to  layer  5,  and  the  middepth  value  assigned  to  layer  3.  The 
layers  were  numbered  from  the  surface  to  the  bottom. 
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Table  2 

Water  Quality  Constituent  Detection  Limits 

Constituent 

Detection  Limit 

Chlorophyll-a 

4.0  mg  m'® 

CBOD5 

1 .0  g  m'® 

PO4-P 

0.02  g  m'^ 

NH4-N 

0.01  g  m'® 

NO2+NO3-N 

0.01  g  m  ® 

c.  Constituent  layer  values  were  calculated  through  arithmetic  averaging. 
The  averaging  procedure  provided  initial  condition  values  for  tempera¬ 
ture  and  DO  for  each  of  the  five  layers  Averaging  provided  values  for 
the  surface,  middepth,  and  bottom  for  the  remaining  constituents;  values 
for  layers  2  and  4  were  calculated  as  the  arithmetic  average  of  adjacent 
layers  (Table  3). 


Table  3 

Initial  Conditions  (Based  on  Laboratory  Analysis) 


Constituent 


Temperature  ®C 


DO  g  Oj  m‘^ 


Algae  g  C  m  ® 


CBOD5  g  Oj  m'^ 


PO,  g  P  m-® 


NH4  g  N  m-® 


NOj+NOa  g  N  m-® 


5 


14.2 


6.3 


0.16 


.5 


0.13 


0.007 


0.06 


Notes:  C/Chl  =  80.0.  Nondetectable  observations  assigned  values  =  0.5  detection  limit. 


Boundary  Conditions 

Boundary  conditions  included  observed  water  quality  at  the  ocean  boun¬ 
dary,  sediment  oxygen  demand  (SOD)  at  the  bottom,  and  light  exchange 
through  the  surface. 

Solar  radiation  at  tte  water  surface  was  computed  from  meteorological 
information.  The  computed  solar  radiation  values  were  used  in  the  algal 
growth  computations.  Meteorological  information  was  acquired  from  the 
National  Oceanic  and  Atmospheric  Administration’s  Tape  Deck  No.  1440 
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WBAN  Hourly  Surface  Observations  for  the  Long  Beach  Airport  (Station 
No.  5085).  Computations  for  solar  radiation  were  based  on  a  Heat  Exchange 
Program  (Eiker  1977). 


Examination  of  measured  SOD  data  (Table  4)  indicated  that  values  varied 
between  1.31  and  2.63  g  Oj  m'^  day'*.  Because  of  the  uncertainties  associated 
with  SOD  (e.g.  measurement  techniques  are  not  well  established  and  SOD 
varies  with  time  and  overlying  DO  concentration),  a  constant,  “representative” 
value  of  2.0  g  O2  m'^  day'*  was  used. 


Table  4 

Sediment  Oxygen  Demand  Measurements 

Station 

Number 

SOD  Measurement 
g  Oj  m‘*  day’’ 

Average  SOD 
g  Oj  m‘*  day  ’ 

Station  Location 

Trial 

2.08 

2.08 

West  LA  Outer 

Harbor 

1-4 

2.63 

1.63 

2.14 

East  Basin  Channel 

1-10 

1.82 

1.82 

Cerritos  Channel 

CO 

1 

1.31 

1.31 

LA  Outer  Harbor 

1-11 

1.63 

1.53 

1.58 

Middle  Harbor 

1-8 

1.75 

1.67 

1.71 

LB  Outer  Harbor 

The  constituents  phytoplankton,  PO4-P,  NH4-N,  NO3-N,  and  CBOD5  were 
set  constant  (Table  5).  The  constituent  values  were  an  average  of  all  boundary 
measurements  (Tekmarine  1987).  Nondetectable  measurements  were  set  equal 
to  one-half  their  detection  limit. 


Table  5 

Temporally  and  Spatially  Constant  Boundary  Conditions 

Constituent 

Value 

Algae  g  C  m'® 

0.2256 

PO4  g  P  m-® 

0.0726 

NH4  g  N  m‘® 

0.0121 

NO3  g  N  m-® 

0.0297 

CBOD5 

0.6300 

The  dissolved  oxygen  and  temperature  at  the  ocean  boundary  were  speci¬ 
fied  by  linearly  interpolating  in  both  space  and  time  the  measured  data. 

Tables  6  and  7  summarize  the  dissolved  oxygen  and  temperature  data  used  for 
the  ocean  boundary. 
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Table  7 

Temperature  (°C)  for  Boundary  Stations  B1,  B2,  and  B3 

Julian  Day 

Layer 

1 

2 

3 

4 

5 

Station  B1 

216.5 

16.1 

15.1 

14.2 

13.5 

12.7 

223.5 

20.4 

19.6 

18.5 

17.8 

17.2 

230.5 

20.4 

20.2 

18.8 

19.0 

15.5 

237.5 

20.7 

20.0 

18.7 

17.9 

16.5 

Station  B2 

216.5 

19.8 

18.4 

16.7 

15.4 

14.4 

223.5 

21.2 

21.4 

19.0 

19.5 

17.6 

230.5 

20.8 

21.8 

19.3 

20.2 

16.6 

237.5 

19.5 

19.0 

18.2 

16.9 

15.8 

Station  B3 

216.5 

21.8 

21.4 

18.2 

18.2 

17.0 

223.5 

22.1 

23.2 

19.6 

21.2 

17.9 

230.5 

21.2 

23.5 

19.8 

21.3 

17.7 

237.5 

20.8 

22.5 

18.8 

20.0 

17.3 
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3  Existing  and  Plan  Flushing 
and  Water  Quality 
Comparison 


The  flushing  studies  consisted  of  insertion  of  a  conservative  tracer  and 
noting  the  movement  and  dilution  of  the  tracer.  The  flushing  studies  provided 
a  qualitative  comparison  between  Existing  and  Pier  400  implementation  condi¬ 
tions.  A  decrease  in  the  flushing  rate  prolongs  the  period  of  time  that  oxygen¬ 
demanding  substances  exert  their  influence  on  the  DO  concentration.  The 
water  quality  simulations  included  the  state  variables  temperature,  phytoplank¬ 
ton,  PO4-P,  NH4-N,  NO3+NO2-N,  CBOD5,  and  DO.  Emphasis  was  on  the  DO. 


Flushing  Comparison 

The  existing  WQM  grid  was  modified  so  that  the  WQM  cells  east  of  the 
access  causeway  coincided  with  the  two  stages  of  the  Pier  400  implementation 
(Plates  13  and  14).  The  modified  existing  WQM  grid  are  referenced  as 
“Existx.”  The  grid  was  modified  because  emphasis  was  directed  toward 
assessment  of  flushing  in  the  lagoon  area  east  of  the  causeway.  A  comparison 
of  Plate  6  (Existing  WQM  Grid  Overlaid)  with  Plate  14  (Modified  Existing 
WQM  Grid  Overlaid)  indicates  that  the  lagoon  area  was  divided  into  two 
WQM  cells  corresponding  to  the  ceU  structure  of  POLA  Stage  1  (Plate  9)  and 
POLA  Stage  2  (Plate  12).  Bold  arrows  were  drawn  on  Plate  13  to  emphasize 
the  location  of  the  two  WQM  cells. 


Flushing  studies 

Two  flushing  studies  were  conducted;  insertion  of  tracer  in  all  WQM  cells 
interior  to  the  breakwater,  and  insertion  of  tracer  in  the  WQM  cell  east  of  the 
proposed  access  causeway  (WQ  Station  X-11). 

The  first  flushing  study  identified  patterns  of  flushing  and  compared  differ¬ 
ences  between  existing  and  Pier  4{X)  implementation  conditions.  The  flushing 
study  with  insertion  of  tracer  in  the  WQM  cell  east  of  the  proposed  access 
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causeway  was  conducted  to  document  changes  in  flushing  due  to  construction 
of  the  access  causeway. 

The  conservative  tracer  was  uniformly  inserted  into  the  five  layers  of  the 
water  column  at  a  concentration  of  10.0  umts  m  ^  at  a  time  corresponding  to 
0000  hr  August  1,  1987.  Tracer  transport  was  simulated  for  30  days.  For  the 
first  flushing  study  (insertion  of  tracer  in  all  the  WQM  cells),  the  boundaries 
were  specified  exterior  to  the  breakwaters.  Tracer  could  exit  the  outer  harbor 
through  the  breakwater  openings,  but  only  water  without  tracer  material  could 
enter  the  outer  harbor  from  the  ocean  boundary. 

Shade  plots  and  graphs  showing  dilution  of  the  tracer  are  displayed  in 
Appendix  A. 


Flushing  results 

Shade  plots  of  the  first  flushing  experiment  are  displayed  in  Plates  A1-A9. 

A  blank  cell  indicates  that  the  simulated  tracer  concentration  was  less  than 
1.0  unit  m  l  Examination  of  Plates  A1-A3  reveals  that  the  pattern  and  magni¬ 
tude  of  tracer  dilution  were  nearly  equivalent  after  5  days  for  Existing,  POLA 
Stage  1,  and  POLA  Stage  2  conditions.  However,  examination  of  shade  plots 
after  15  days  (Plates  A4-A6)  and  25  days  (Plates  A7-A9)  reveals  that 

a.  POLA  Stage  1  and  Stage  2  conditions  inhibit  flushing  in  the  LA  Outer 
Harbor,  Fish  Harbor,  Seaplane  Lagoon,  and  Main  Channel. 

b.  Flushing  patterns  are  similar  for  both  POLA  Stage  1  and  Stage  2. 

It  is  important  to  note  that  the  insertion  of  tracer  in  all  WQM  cells  provides 
information  on  global  flushing  but  provides  little  information  on  flushing  at 
specific  sites.  For  example,  the  shade  plots  indicates  that  high  concentration  of 
tracer  remains  in  the  Seaplane  Lagoon  throughout  the  30-day  simulation 
period.  Similarly,  the  shade  plots  reveal  that  high  concentrations  of  tracer 
remain  in  water  adjacent  to  the  Seaplane  Lagoon.  Any  tracer  flushed  from  the 
Lagoon  would  be  replaced  with  adjacent  water  with  high  tracer  concentration. 

Results  of  the  second  flushing  study,  insertion  of  tracer  in  the  WQM  cell 
east  of  the  proposed  causeway,  is  displayed  using  a  time  series  plot  showing 
tracer  concentration  with  time  (Plate  AlO).  Examination  of  Plate  AlO  indi¬ 
cates  that  flushing  is  equivalent  in  both  POLA  Stage  1  and  POLA  Stage  2 
conditions.  Existing  conditions  required  2.5  days  for  the  tracer  to  be  diluted 
90  percent,  while  POLA  Stage  1  and  2  conditions  required  7.5  days.  Exami¬ 
nation  of  the  WQM  grid  for  the  existing  condition  (Plates  13  and  14)  shows 
that  advection  to  cells  west  of  the  insertion  point  was  possible  because  of  the 
absence  of  the  planed  access  causeway.  A  shade  plot  (Plate  All)  of  a  snap¬ 
shot  taken  2.5  days  after  tracer  insertion  in  the  existing  condition  shows  that 
tracer  was  transported  to  cells  west  of  the  insertion  point. 
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The  access  causeway  prevented  tracer  advection  to  cells  west  of  the  inser¬ 
tion  point.  Examination  of  the  shade  plot  for  the  POLA  Stage  2  condition 
(Plate  A 12)  reveals  no  tracer  west  of  the  access  causeway. 


Water  Quality  Comparison 

Two  techniques  were  used  to  compare  measured  data  with  Base,  Existing, 
POLA  Stage  1,  and  POLA  Stage  2  simulations:  time  series  plots,  and  vertical 
profile  plots. 

Snapshots  of  the  simulations  were  made  to  correspond  to  the  time  of  collec¬ 
tion  of  measured  data.  The  snapshots  were  then  displayed  as  time  series  and 
vertical  profile  plots.  The  location  and  naming  convention  of  points  used  for 
comparison  (displayed  in  Plates  15-18)  is  described  below. 


Plate 

Plot 

15 

Water  Quality  Stations  -  Base 

16 

Water  Quality  Stations  -  Existing 

17 

Water  Quality  Stations  -  POLA  Stage  1 

18 

Water  Quality  Stations  -  POLA  Stage  2 

The  prefix  “I”  refers  to  the  interior  stations,  and  the  prefix  “B”  refers  to  the 
boundary  stations  sampled  by  Tekmarine.  The  prefix  “X”  represents  the  extra 
locations  selected  for  more  complete  coverage  of  the  study  area.  The  stations 
in  each  of  the  simulations  correspond  to  the  same  physical  location;  Plates  IS¬ 
IS  provide  guidance  in  locating  station  locations  relative  to  plan  morphology. 
The  water  quality  stations  and  the  centroid  of  the  water  quality  cell  used  to 
represent  the  stations  are  listed  in  Table  8.  The  centroid  is  represented  as 
hydrodynamic  grid  cell  coordinates  (i,j). 


Time  series 

The  water  quality  constituents  phytoplankton,  PO4-P,  NH4-N,  NO3-1-NO2-N, 
CBOD5,  and  DO  are  displayed  as  time  series  for  the  surface  (layer  1)  and 
bottom  (layer  5)  layers.  Snapshots  were  made  at  1200  hr  August  4  (Julian 
Day  216),  August  11  (Julian  Day  224),  August  18  (Julian  Day  230),  and 
August  25  (Julian  Day  237).  The  snapshot  values  were  linearly  interpolated 
between  snapshot  dates  and  displayed  in  the  time  series  plots  in  Appendix  B. 
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Table  8 

Water  Quality  Stations 

Station 

Number 

HM  Grid 

Description 

1 

J 

I-01 

14 

48 

Main  Channel 

1-02 

26 

48 

TITP  Outfall 

1-03 

26 

37 

LA  Outer  Harbor 

1-04 

26 

27 

East  Basin  Channel 

1-05 

49 

60 

Middle  Harbor 

1-06 

57 

69 

Back  Channel 

1-07 

82 

62 

LA  River  Mouth 

1-08 

75 

29 

LB  Outer  Harbor 

1-09 

98 

20 

SE  of  Island  Freeman 

1-10 

38 

81 

Cerritos  Channel 

1-11 

57 

60 

Middle  Harbor 

B-01 

26 

18 

Angels  Gate  Entrance 

B-02 

73 

18 

Queens  Gate  Entrance 

B-03 

108 

13 

E  LB  Breakwater 

X-01 

22 

90 

LA  West  Basin 

X-02 

40 

91 

Consolidated  Slip 

X-03 

12 

80 

LA  Southwest  Slip 

X-04 

64 

82 

LB  Channel  Two 

X-05 

64 

77 

LB  Channel  Three 

X-06 

31 

64 

Seaplane  Lagoon 

X-07 

38 

64 

Naval  Basin 

X-08 

6 

44 

LA  West  Channel 

X-09 

68 

63 

LB  East  Basin 

X-10 

44 

45 

LB  Outer  Harbor 

X-11 

35 

64 

Access  Causeway 

X-12 

3 

25 

San  Pedro  Breakwater 

Vertical  profiles 

Examination  of  the  time  series  plots  revealed  DO  differences  between 
measured  data,  Existing,  POLA  Stage  1,  and  POLA  Stage  2  conditions  for 
Stations  I-Ol,  1-02, 1-07,  and  X-12.  The  water  quality  snapshots  for  these 
stations  were  displayed  using  vertical  profiles.  Station  X- 11,  the  location  east 
of  access  causeway,  was  included  in  the  vertical  profiles.  The  vertical  profile 
plots  are  displayed  in  Appendix  C. 
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Results 


Examination  of  water  quality  simulation  results  of  Base  and  Existing  condi¬ 
tions  indicates  that  Base  and  Existing  conditions  were  nearly  equivalent. 

The  measured  and  simulated  results  for  the  Base  and  Existing  conditions 
were  in  general  agreement  except  for  Stations  1-2  (Plates  B2,  B28,  Cl,  C6, 
Cll,  and  C16)  and  1-7  (Plates  B7,  B33,  C3,  C8,  C13,  and  C18).  Station  1-2 
was  located  at  the  Terminal  Island  Treatment  Plant  (TITP)  outfall,  and  the 
measured  constituent  values  reflected  the  effluent  conditions  during 
August  1987.  All  the  model  simulations  did  not  include  the  TITP  effluents. 
The  TITP  outfall  effluent  was  not  modeled  because  other  organizations  respon¬ 
sible  for  the  TITP  outfaU  have  modeling  responsibility  for  impacts  separate 
from  Pier  4(X).  The  TITP  outfaU  wiU  be  relocated  in  the  future  regardless  of 
harbor  expansion  plans. 

Examination  of  the  plots  for  Station  1-7  reveals  that  the  Los  Angeles  River 
is  contributing  some  flows  to  the  bay  that  were  not  modeled.  The  water 
quality  sampling  program  conducted  by  Tekmarine  revealed  less  saline,  as  well 
as  nutrient-,  CBOD5-,  and  chlorophyU-enriched  water  at  Station  1-7.  The  pro¬ 
file  plots  for  Station  1-7  (Plates  C3,  C8,  C13,  and  C18)  indicate  that  the  Los 
Angeles  River  may  have  deposited  organic-enriched  sediments  at  the  location, 
resulting  in  increased  sediment  oxygen  demand.  An  increased  SOD  would  be 
indicated  by  DO  gradients  with  maximum  DO  at  the  surface  decreasing  to 
minimum  DO  at  the  bottom.  Measured  DO  with  surface  concentrations 
exceeding  10.0  and  bottom  concentrations  less  than  4.0  g  DO  m'^  were  mea¬ 
sured  at  Station  1-7. 

A  comparison  of  Existing,  POLA  Stage  1,  and  POLA  Stage  2  water  quality 
time  series  and  profile  plots  reveals  that  maximum  differences  in  DO  occurred 
at  Station  I-Ol  (Main  Channel),  1-02  (TITP  Outfall),  and  X-12  (San  Pedro 
Breakwater).  An  examination  of  the  time  series  plots  indicates: 

a.  Maximum  differences  occurred  the  second  week  of  simulation  (Julian 
Day  223.5). 

b.  Simulated  DO  was  nearly  equivalent  at  the  beginning  and  end  of  the 
simulation  period.  The  similarity  at  the  beginning  of  the  simulation 
period  was  due  to  the  same  initial  conditions.  The  simulated  differ¬ 
ences  the  second  week  may  be  a  transient  response,  and  the  similarity 
at  the  end  of  the  simulation  period  may  reflect  that  water  quality 
conditions  following  Pier  400  implementation  are  similar  to  existing 
conditions. 

c.  Maximum  differences  between  Existing  and  the  Pier  400  implementa¬ 
tions  in  DO  occurred  in  the  bottom  layer  at  Station  1-02  at  Julian  Day 
237.5.  The  simulated  existing  DO  was  8.6  while  POLA  Stage  1  was 
7.2  and  POLA  Stage  2  was  7.1  g  DO  m'^.  The  simulated  bottom  layer 
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DO  differences  between  Existing  and  the  Pier  400  implementations 
were  near  1.5  g  DO  m'^;  in  contrast,  surface  layer  differences  were  less 
than  1.0  g  DO  m"^. 

d.  The  maximum  simulated  differences  in  DO  at  Station  I-Ol  occurred  at 
the  surface.  Existing  surface  DO  was  8.6  g  DO  m'^  while  POLA  Stage 
1  was  8.1  and  POLA  Stage  2  was  7.9  g  DO  m'^.  The  simulated  differ¬ 
ences  were  less  than  1.0  g  DO  m'^. 

e.  The  Pier  400  implementations  decreased  the  simulated  surface  and  bot¬ 
tom  DO  at  Station  X-12.  Surface  DO  decreased  about  1.1,  and  bottom 
DO  decreased  0.5  g  DO  m'^. 

/.  Examination  of  the  simulated  DO  at  Station  X-11,  located  at  the 
embayment  east  of  the  access  causeway,  revealed  that  simulated  DO 
was  greater  in  the  Pier  400  implementations  than  in  the  Existing  condi¬ 
tions  (Plates  B25,  B51,  C4,  C9,  C14,  and  C19).  Maximum  differences 
occurred  the  second  week  of  simulation,  with  both  surface  and  bottom 
DO  in  the  Pier  400  implementations  exceeding  existing  conditions  by 
0.3  g  DO  m  l 

Generally,  the  simulated  DO  differences  between  the  Existing  and  Pier  400 
implementation  conditions  were  minimal.  Simulated  DO  decreased  slightly  in 
the  Los  Angeles  Outer  Harbor  and  Main  Channel  following  Pier  400  stage 
implementations.  The  magnitude  of  the  DO  differences  increased  with  ftirther 
plan  implementation;  that  is,  POLA  Stage  2  resulted  in  lower  DO  concentra¬ 
tions  than  POLA  Stage  1,  and  POLA  Stage  1  resulted  in  lower  DO  than  Exist¬ 
ing.  The  only  model  input  differences  between  Existing  and  the  Pier  400 
implementations  were  changes  in  basin  morphology  due  to  dredge  and  fill. 
Simulated  differences  were  due  to  changes  in  circulation  and  flushing  as  indi¬ 
cated,  with  inhibited  flushing  in  the  LA  Outer  Harbor,  Fish  Harbor,  Seaplane 
Lagoon,  and  Main  Channel  following  Pier  400  implementation.  The  flushing 
studies  indicated  a  change  in  circulation  and  flushing. 


Chapter  3  Existing  and  Plan  Flushing  and  Water  Quality  Comparison 


15 


4  Summary  and  Conclusions 


A  numerical  water  quality  model  was  applied  to  San  Pedro  Bay,  and  simu¬ 
lation  results  for  Existing  and  POLA  Stage  1  and  POLA  Stage  2  implementa¬ 
tions  of  the  Pier  400  port  facility  were  compared.  The  focus  of  the  study  was 
to  compare  flushing  and  DO  resources. 

The  flushing  studies  consisted  of  insertion  of  a  conservative  tracer  and 
noting  movement  and  dUution  of  the  tracer.  A  flushing  study  with  insertion  of 
tracer  in  all  WQM  ceUs  indicated  that  the  Pier  400  implementation  inhibited 
global  flushing  on  the  west  side  of  the  Los  Angeles  Outer  Harbor. 

The  insertion  of  tracer  in  the  WQM  ceU  east  of  the  proposed  access  cause¬ 
way  revealed  that  construction  of  the  causeway  resulted  in  a  300  percent 
decrease  in  the  rate  of  dilution. 

The  water  quality  results  indicate  that  the  simulated  DO  is  nearly  equivalent 
between  Existing  and  the  Pier  400  conditions.  Maximum  differences  in  DO 
occurred  in  the  west  side  of  the  Los  Angeles  Outer  Harbor.  Decreases  in  DO 
were  noted  in  the  Main  Channel,  the  prior  TITP  effluent  location,  and  in  the 
Outer  Harbor  interior  to  the  San  Pedro  Breakwater.  The  maximum  decrease  in 
DO  occurred  at  the  prior  TITP  effluent  location,  located  north  of  the  Pier  400 
landfill.  The  simulated  bottom  layer  DO  for  POLA  Stage  2  was  7.1  g  DO  m' , 
a  decrease  of  nearly  1.5  g  DO  m‘^  from  Existing  conditions.  The  simulated 
DO  concentration  of  7.1  g  DO  m'^  is  near  90  percent  of  DO  saturation  at  the 
measured  temperature  and  salinity.  The  simulated  DO  concentration  is  ade¬ 
quate  to  maintain  the  existing  aquatic  biota.  Surface  layer  differences  were  less 
than  1.0  g  DO  m‘^. 

Examination  of  the  simulated  DO  in  the  embayment  east  of  the  access 
causeway  revealed  that  simulated  DO  was  greater  in  the  Pier  400  implementa¬ 
tions  than  in  the  Existing  conditions.  DO  in  the  Pier  400  implementations 
exceeded  Existing  conditions  a  maximum  of  0.3  g  DO  m'^. 

Examination  of  results  for  Existing,  POLA  Stage  1,  and  POLA  Stage  2 
indicates  that  simulated  DO  exceeds  6.0  g  DO  m'^  for  all  sampled  stations. 
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